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We consider the non- equilibrium ferromagnetic transition of a mesoscopic 
sample of a resistive Stoner ferromagnet coupled to two paramagnetic leads. 
The transition is controlled by either the lead temperature T or the transport 
voltage V applied between the leads. We calculate the T and V dependence 
of the magnetization. For systems with a flat density of states we find within 
mean-field theory that even at finite bias the magnetization does not depend 
on the position along the sample axis, although the charge density and other 
quantities do vary. This may be relevant for possible spintronics applications. 
In addition, we establish a generalized control parameter in terms of T and 
V which allows for a universal description of the temperature- and voltage- 
driven transition. 
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1. INTRODUCTION 

Phase transitions in interacting electron systems far from thermody- 
namic equilibrium are an exciting, relatively young field of research, both 
for fundamental theoretical reasons and for application purposes. While the 
theory of equilibrium phase transitions is well developed, a theory for phase 
transitions out of equilibrium is much more difficult to formulate, largely 
because in the most general case the thermodynamic states of the system 
are not characterized by a minimum of the free energy and because often 
the temperature T as a state variable driving the transition is not even 
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defined. On the other hand, nonequilibrium phase transitions are of high 
technological potential, since a mesoscopic electron system can be driven 
out of equilibrium in a highly controlled way, and controlling a nonequilib- 
rium phase transition by externally applied fields can be envisaged. It opens 
up the prospect of fast switching between different magnetic or conducting 
states of a mesoscopic system, e.g. by applying a transport bias voltage or 
by an applied laser field, with possible spintronics applications. 

A thermodynamic phase transition in equilibrium is in general controlled 
by the minimum of the free energy 

T = U-TS, (1) 

i.e. by the balance of the internal energy U and the entropy S. The depen- 
dence on the control parameter temperature favors an ordered phase at low 
T and a disordered phase at high T. In a stationary nonequilibrium situa- 
tion internal energy and entropy are statistically still well-defined quantities. 
One may, therefore, conjecture that a relation analogous to Eq.[T] holds even 
in this case. The different phases of a stationary non-equilibrium system 
would then again be determined by the minimum of the free energy. How- 
ever, in the above relation (TjQ) the entropy coefficient T would have to be 
generalized, since temperature T is in general not defined in nonequilibrium. 
Based on general thermodynamic arguments we propose in the present ar- 
ticle a definition of such a generalized control parameter r which allows to 
describe electron systems driven out of equilibrium by an applied stationary 
bias voltage V in a unified way. 

To be specific, we consider the voltage-driven magnetic phase transition 
of a Stoner ferromagnet. Such a system can be realized as a mesoscopic 
sample of a magnetic metal, contacted by two paramagnetic electron reser- 
voirs or leads at chemical potentials (jll and fin, respectively, as seen in Fig. 
[Th). It has recently been demonstrated that the fabrication of such devices, 
which sustain a high nonequilibrium bias, is technologically feasible using 
metallic nanobridges with high lead-to-bridge aspect raticP or using resis- 
tive nanowiresP As will be seen below, the ferromagnetic transition in the 
mesoscopic sample can be controlled both by the reservoir temperature T 
and by the voltage V = fii — fJ-R applied between the reservoirs. The mag- 
netic state of the sample may be probed by a third, ferromagnetic electrode, 
attached to the sample via a tunnel junction. This electrode stays in equilib- 
rium with itself, independent of the transport voltage V, so that it remains in 
the ferromagnetically ordered state throughout the transition of the sample, 
and its tunneling conductance depends strongly on the magnetization in the 
sample. A related proposal of controlling magnetism by a nonequilibrium 
transport bias has recently been put forward in Ref. [3]. However, in that 
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proposal the magnetization was induced in a paramagnetic metal sample 
by ferromagnetism in the leads, and the sample was assumed to be zero- 
dimensional in the sense that it was smaller than the electrons' elastic mean 
free path £, implying ballistic transport and a position independent electron 
distribution function in the metal sample. In the present work we consider 
the more realistic case of a resistive, ferromagnetic sample whose lateral size 
L is greater than the elastic mean free path, L > £, but still smaller than the 
inelastic and spin relaxation lengths, L < £i n ,£ s . In this case the transport 
is diffusive, leading to a nonequilibrium electron distribution function which 
depends on the location along the sample. We develop a mean field theory 
for the temperature- or voltage-driven Stoner transition. 

One might think that this transition, driven by the transport voltage 
at reservoir temperature T = might be a quantum phase transition, since 
it is controlled by a parameter different from the temperature. However, 
at any finite voltage V ^ the nonequilibrium quasiparticle distribution 
is such that there exists a finite phase space for inelastic scattering (see, 
e.g., Section 2.1), so that the quantum dephasing rate in the sample, 7^, 
remains finite even for reservoir temperature T = 0. Thus, critical diver- 
gencies of correlation functions at the phase transition are dominated by 
classical, statistical fluctuations rather than quantum fluctuations, i.e. the 
transition belongs to the universality class of classical rather than quantum 
phase transitions, in contrast to the assertion of Ref. [4]. Remarkably, the 
authors of Ref. [4] do, however, find that the transition is governed by the 
classical critical exponents. In fact, as will be seen in Section 3, the roles of 
V and T are completely analogous, as far as thermodynamic properties are 
concerned. 

This paper is organized as follows. In Section 2 we define the Stoner 
model in a steady-state nonequilibrium situation with finite bias. The quan- 
tum Boltzmann equation for a resistive wire with a Stoner mean-field in- 
teraction is derived in Section 2.1, and the Stoner mean- field equations in 
the presence of finite bias are presented in Section 2.2, along with their 
numerical solution for various temperatures and voltages. In Section 3 we 
provide general thermodynamic arguments which allow us to introduce a 
control parameter t(T,V), depending on the reservoir temperature T and 
the bias voltage V, as a generalization of the temperature to the voltage- 
biased non-equilibrium situation, as far as thermodynamic properties are 
concerned. We demonstrate explicitly that this provides a universal descrip- 
tion of the voltage-driven Stoner transition. Some concluding remarks are 
given in Section 4. 
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Fig. 1. a) Schematic view of a magnetic transitor. The ferromagnetic meso- 
scopic sample is shown in light grey. It is attached to paramagnetic leads 
kept at chemical potentials and hr, respectively. A magneto-tunneling 
current depending on the sample magnetization can be extracted through 
the third electrode FM. The position coordinate along the sample of length 
L is denoted by s, as shown in the figure, b) The Keldysh Hartee-Fock di- 
agrams defining the Stoner mean-field theory in stationary nonequilibrium. 
Solid lines denote the electron propagators, wavy lines the screened Coulomb 
interaction U. 

2. VOLTAGE-DRIVEN STONER FERROMAGNETISM 

As a simple model for a ferromagnetic transition in an itinerant, diffusive 
electron system we consider electrons on a lattice with dispersion e^, with 
a random potential Vj induced by static impurities on lattice sites j and a 
strongly screened, local Coulomb repulsion U between electrons on the same 
site, 

h = £ ekcL^T + E v A c * + U Y1 c ^ c jAi c n ■ ( 2 ) 

kcr j cr j 

Here , c k(j denote the electron creation and destruction operators with 
momentum k and spin a = ±1/2, and cj CT = ^ k exp(zkxj)c kcr is the creation 
operator at site j. We assume that any antiferromagnetic instabilities that 
could be induced by nearest neighbor hopping on a cubic lattice are frus- 
trated, e.g., by an appropriate lattice structure or by a next-nearest neighbor 
hopping included in the band structure e^. 

2.1. Nonequilibrium Distribution in a Diffusive Ferromagnet 

Any theory of an interacting system out of equilibrium requires knowl- 
edge of the quasiparticle distribution faxi^) as a function of the particle 
energy to. In a paramagnetic, diffusive system at bias V, if inelastic relax- 
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ations can be neglected, f ax {^) is knowiJSKEE i Q have a double-step form 
according to the two Fermi edges in the reservoirs, and to depend linearly on 
the position x = s/L along the sample axis (see Fig. [TK))- In a ferromagnet, 
however, the Stoner splitting may in general depend on x as well, inducing 
nontrivial position and spin dependences to f ax {w) even in the absence of 
inelastic or dynamical spin flip processes. In the following we derive the dis- 
tribution function for the case of a mesoscopic ferromagnet, where ineleastic 
and spin flip processes can be neglected. 

The kinetic equation (quantum Boltzmann equation) for f a x(^) is ob- 
tained as the off-diagonal part of the Keldysh Dyson equation for the lesser 
Green's function G^(u, x, x')P After transforming to center-off-mass and 
relative coordinates, r = (x + x') /2, Ax = x — x', taking Ax to be confined 
in a small volume compared to the scale on which the external fields vary, 
and Fourier transforming to momentum space with respect to Ax within 
that volume, the resulting quantum Boltzmann equation acquires the form 
of a continuity equation in phase space, 



^ + v k • V r + yE(r) • V fc 



G< r (k,u)=l{f ar (k)} 



(3) 



where G< r (k, uj) = 27ri/ OT (k).A OT (k, to) is expressed in terms of the distribu- 
tion /o-r(k) and the spectral density A ar (k,uj) of momentum k at position 
r in the sample, = dr/dt = de^/dk is the electron group velocity and 
— eE(r) = hdk/dt the force acting on the electrons due to an external field 
E(r). The latter includes both, the external bias voltage and the electrostatic 
potential generated by the random impurities. Since we consider a station- 
ary situation and assume that the sample is smaller than the inelastic and 
the spin relaxation lengths, the explicit time derivative and the collision in- 
tegral I {f ar (k)} vanish in Eq. 0. After summation over k the force term 
vanishes as well. This is because of the relation 

£ V k • E(r)/ r (k)A,r(k, w) = <f E(r)/„(k)A„(k, uj) ■ dS k = . (4) 



kcr 



Here the second term implies an integration over the surface of the 1st 
Brillouin zone, where / CTr (k) is either or 1, i.e. f a r{— k) = far (k) and 
dS_k = — dS^. The same k-summation also introduces the particle and the 
current density in Eq. ([3]), 



p ar (uj) = ^/^(kjArrfk,^) 
k 

JcrrH = ^2 V k f ar (k)A ar (k,Uj) , 



(5) 
(6) 
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which in a diffusive system with diffusion coefficient D are related by Fick's 
law, 

j CTr M = -DV r p r (u) . (7) 

In this way the diffusive Boltzmann equation in a ferromagnet takes the form 

-DV*[f„(u)N ar (u)] = 0, (8) 

where N ar (uj) = ^ k A ar (k, u) is the Stoner split electronic density of states 
(DOS) at position r in the sample. Note that in a paramagnetic system 
No-r^) is in general position independent and can be dropped from Eq. 
(JH). In a ferromagnet, however, the spectral function contains the Stoner 
selfenergy S,j r (see Section 2.2), which depends in general on the position 
r = s, 

A ax (k,uj) = Im— , (9) 

where s is the coordinate along the sample axis (see FigQJ-)) and x = s/L. 
Taking into account only the linear dependence of the DOS on x, with 
N' a (u>) = dN ax (u)/dx\ x=0 .5 = -[dN ax {uj)/duj ■ dT, ax /dx]\ x=0 , 5 and N a (tn) = 
N ax {uj)\ x= Q^, the quantum Boltzmann equation in a resistive ferromagnet 
takes the final form 

dx* lU N a (u) dx ~ U ' UUj 

The second term in this equation implies a selfconsistent coupling to the 
Stoner mean-field equations via the spin and position dependent DOS. How- 
ever, Eq. (fTU|) can readily be solved analytically, once N^(uj) and N a (u>) are 
determined. This will be done in a forthcoming publicationP In the present 
paper we wish to focus on universality related issues, which can be discussed 
as well for the particle-hole symmetric case. Therefore, we will assume in 
the following that the DOS is flat within the energy window defined by the 
applied bias voltage V. In this case N' a (uj) = 0, and the kinetic equation 
(jlOP reduces to the form familiar from noninteracting diffusive systems. In 
particular, f ax (u) is then spin independent, since the boundary conditions 
on f ax (oo) are spin symmetric, i.e. both spin species in the sample couple to 
the same paramagnetic reservoirs. 

faxH = f- ax {u) = f x (uj) for N' a (uj) = . (11) 



The Boltzmann equation fjlOj) is easily solved for this case with the boundary 
conditions that in the left (right) lead there is an equilibrium distribution 
with chemical potenial hl = +eV/2 and fj,R = —eV/2, respectively, 

fx(u) = xf^(co + ^) + (1 - z)/ (0) (^ - ^0 , (12) 
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where f(°\iv) = l/[exp(u)/kBT) + 1] is the Fermi distribution and T the 
reservoir temperature. 

2.2. Stoner Mean-Field Equations in Stationary Nonequilibrium 

For simplicity we assume an elliptic DOS (per spin) of the disordered 
system without interactions, Nq(u> + //q) = (2/-7r)y 1 — {u + /^o) 2 > where the 
single-particle energy uj is measured relative to the equilibrium chemical 
potential [iq, and all energies are given in units of the half bandwidth B. 
To realize a flat DOS for energies within the voltage window, —eV/2 < lo < 
eV/2, fiQ is placed at the band center, fiQ = 0, and eV <C B. As discussed 
in Section 2.1, the distribution function is then decoupled from the solution 
of the Stoner model and is given by Eq. f)12[> . 

The Hartree-Fock selfenergy for this model is represented by the 
Keldysh diagrams shown in Fig. QJ)). Note that all vertex points of the 
instantaneous interaction U are on the same branch of the Keldsyh contour, 
as indicated by the "— " signs, and that the equal-time operator products 
are normal-ordered. The interal equal-time propagators are, therefore, the 
lesser Green's functions G < (a;)PThe Hartree-Fock selfenergy then reads, 

E a (T,V,x) = U(n- a ) 

= U [ duf x {u)N [u + iio-V-*(T,V,x)] , (13) 



where (n a ) is the conduction electron number per site with spin a. At 
T = the integral can be evaluated, and the V and x dependence can be 
made explicit, 



U 

TT 



x ^J)_ ) _ CT ^1 — (il-^o-) 2 + arcsin(r2_ ) _ cr )^ 



(14) 



+ (l-x)[ n+ a Jl - (tt+-a) 2 + arcsin(0 + _ CT ) + vr 



with n± )(T = + hq - S ff (0, V, x). The selfenergy can be separated in a 
spin independent and a spin dependent part, T, a = So + crAS, with 

AS = St-S^ (15) 

So = ^(S T + S i ) = |(n) a; . (16) 

It is now crucial that the definition (I12p of the nonequilibrium distribution 
fixes the zero of the energy to be in the middle between the two Fermi 
steps for each position x along the sample. Therefore, in Eqs. (|13p . (j!4[) 
the interaction-induced, spin independent shift of the energy Sq must be 
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Fig. 2. The magnetization per site is shown as a function of T for various 
values of the transport voltage V as indicated. V is given in units of the 
equilibrium Curie temperature T c q. The screened Coulomb interaction was 
chosen such that T c0 = 0.01 B (i.e. UN o (0) = 1.0001669). The small dips in 
the m(T) curves near T = are a consequence of the slightly non-constant 
DOS for high bias V. 

absorbed into the (equilibrium) chemical potential for each x: fiQ— So — > Ho- 
The selfenergy then consists only of the spin dependent part, E CT = erAE. 
Eq. (fT5|) . together with Eqs. {JT2]) and (fT3j) constitute the unconstrained, 
selfconstistent Hartree-Fock equations which determine the magnetization 
per site, m = g^ueAE/i/, g and being the Lande factor and the Bohr 
magneton, respectively. We note that for the present case of a flat DOS it 
follows immediately from Eqs. (I13p . (|15fl . (I16p and from the linearity of f x (oj) 
with respect to x, Eq. (|12p . that (i) AS, and hence the magnetization, are 
independent of the position x and that (ii) the total particle density (n) x 
is a linear function of x. The latter is consistent with Fick's law (|7|) and 
the fact that the stationary current j x in the sample is homegenous. The 
x independence of the magnetization m remarkably implies that the total 
sample magnetization can be switched as a whole by a transport voltage V, 
although at finite bias the electron density (n) x and other quantities vary 
along the sample. This may be of technological relevance for spintronics 
devices, although this property does not persist exactly for non-constant 
DOSP The mean-field solutions for the magnetization m as a function of 
reservior temperature T are shown for various applied transport voltages V 
in Fig. [2] and as a function of V for various T in Fig. [3j For these data 



. 


1 1 1 1 1 1 

V = 0T c0 

0.5 


- 


?~0 — 

2.5^ — 


^~~-\^^^S^ l.O 




_ 31) -^.^^ 
3^2~~~~~~-^^ 












^3^6^ ( ^ 







Universality in Nonequilibrium Phase Transitions 



0.02 




0.015 - 



w 0.01 - 



0.4 



B 



0.5 



0.005 - 



0, 







1 



2 



3 



V/T c0 



Fig. 3. The magnetization per site is shown as a function of V for various 
values of T, as indicated. The parameter values used are as in Fig. El 

the strength of the screened Coulomb interaction U has been chosen such 
that the Stoner criterion, UNq(0) > 1, is satisfied and that the equilibrium 
Curie temperature is T c q = 0.01 B, corresponding to T c q 300 K for a half 
bandwidth of B = 3 eV. 



The fact that the Stoner transition can be controlled by both, temper- 
ature T and transport voltage V, the similarity of the magnetization curves 
m(T) and m(V) in Figs. [2] and [3l and the independence of the magneti- 
zation of the position x in the sample for a flat DOS suggest that there 
is a universality in this transition in terms of the two control parameters 
T and V. Indeed, the complete thermodynamics of any system is deter- 
mined, once the internal energy IA or the free energy T (with appropriate 
state variables) is known. At least the internal energy is well-defined also 
in a stationary nonequilibrium situation. For any stationary electron sys- 
tem in or out of equilibrium, IA is uniquely determined by the quasiparticle 
distribution function f x (u>). The latter is, in turn, uniquely defined by the 
applied transport voltage V (and the concomitant boundary conditions), if 
the system is coupled to one or more reservoirs to ensure stationarity, and 
by the reservoir temperature T. This suggests that it should be possible to 
define a single, unified state variable, involving both T and V, in terms of 
the internal energy. 



3. UNIVERSALITY 
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X(T,V)/T c0 

Fig. 4. Collapse of all magnetization data contained in Figs. [2] and [3] onto 
a single curve, depending on t(T, V) (Eq. [TBj) only. The data for different 
parameter values of T and V are displayed here using different symbols. 



Therefore, we put forward the following procedure. The temperature of 
an equilibrium system is usually measured by bringing it in heat contact with 
a reference system whose internal energy is known as a function of T, like 
an ideal gas (ideal gas thermometer). Since the ideal gas is not suitable for 
applying a transport voltage, we propose to use a diffusive, noninteracting 
electron g reference system for any voltage-biased, interacting system, 

like, e.g., a mesoscopic Stoner ferromagnet, Eq. (|2|). As before, we will 
restrict ourselves here to the case that the DOS is constant within the voltage 
window, iVo(O), and defer the general case to a forthcoming paper P The 
change of the internal energy density of this reference system as a function 
of T and V is well-defined and is readily calculated as, 



AU 



U(T, V) — U(0, 0) 



iVo(O) 



yiVo(O) 



(17) 



( k BTf + — 2 {eVf 



Remarkably, AU does not depend on the position x where the energy density 
is considered, as long as the DOS is flat, in accordance with the results of 
Section 2.2. The expression (|17p can now be used to extract a generalized 
state variable, t(T,V), that uniquely determines the internal energy of the 
reference system, and thus the thermodynamics of any stationary electron 
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system controlled by the state variables T and V, 

r(T,V) = ^Jr^+^V 2 , (18) 

where C = Tr 2 kg/3e 2 is the Lorentz number, and AU = (tt 2 /6)Nq(0)t 2 . 
Note that the form of the expression for r(T, V) is equivalent to the effective 
temperature of a voltage-biased, hot electron system that has equilibrated 
due to strong inelastic processes P However, in our line of arguments above no 
assumption has been made about whether or not there are inelastic processes 
involved. This supports the universality of the definition Eq. (|18p . 

We have tested whether the Stoner transition can be described by 
r(T, V) in a universal way. In Fig. O we have plotted all the data con- 
tained in Figs. [2] and [3] as a function of t(T,V). Indeed there is a perfect 
collapse of all data onto a single curve, as expected. 



4. CONCLUSION 

Even the simple Stoner model of a ferromagnetic transition can ex- 
hibit interesting behavior when considered in a non-equilibrium situation. 
We have used this model to define a generalized control parameter r(T, V) 
which allows to describe the voltage and temperature dependence of the 
Stoner transition in a unified way. The model has been evaluated within 
mean-field theory where the critical exponent of the magnetization in de- 
pendence of both T and V is a = 1/2. We emphasize that the generality of 
our thermodynamic arguments that led to the definition of t(T.V) strongly 
suggest that it not only allows for a universal description of the magnetiza- 
tion but also of the complete thermodynamics of the Stoner model, and that 
this holds true even beyond mean field theory. In particular, correlations will 
certainly change the value of a, but it will remain the same for the T and the 
V dependence. Moreover, the thermodynamic arguments suggest that not 
only Stoner ferromagnets but the thermodynamics of any interacting elec- 
tron system controlled by the reservoir temperature T and a nonequilibrium 
transport voltage V can be described by r(T, V) in a universal way. Con- 
versely, it is obvious that dynamical and spectral properties are influenced 
by T and V in quite different ways, which cannot be captured only by the 
T and V dependence of the internal energy. Therefore, no unified behavior 
is expected for these quantities. 

As mentioned above, the magnetization m will acquire an x-dependence, 
as soon as the density of states in the sample as function of energy is not flat. 
In this case the interaction-induced correlation length will also influence the 
spatial dependence of m. The thermodynamics of voltage-driven systems, 



J. Kroha, M. Arnold, and B. Griepernau 



the influence of correlations in general and the position dependence of the 
magnetization in particular will be subjects of further research. 
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